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ABSTRACT 


Sound  waves  propagating  through  the  oceans  are  refracted  by  internal  waves.  In 
the  shallow  waters  of  the  continental  shelves,  an  additional  downward  refraction  of  sound 
waves  due  to  internal  waves  can  cause  them  to  interact  more  often  with  the  seabed, 
resulting  in  additional  energy  from  the  sound  waves  being  dissipated  into  the  seabed. 

This  study  investigates  how  internal  waves  affect  sound  propagation  on  the 
continental  shelves.  It  first  quantified  the  types  of  internal  waves  on  the  continental  shelf 
of  California,  near  Point  Sal,  using  data  collected  from  a  field  experiment.  Next,  the 
effects  that  these  internal  waves  have  on  sound  propagation  were  quantified  via 
simulations  using  a  ray  theory  acoustic  model. 

The  results  showed  that  internal  waves  in  the  experiment  area  were  largely 
generated  by  tidal  forcing.  Compared  to  simulations  without  internal  waves,  simulations 
accounting  for  the  effects  of  internal  waves  resulted  in  higher  sound  energy  loss,  as 
internal  waves  tend  to  cause  sound  waves  to  strike  the  seabed  at  steeper  angles  and  over 
shorter  distances.  Thus,  to  enable  a  more  accurate  assessment  of  underwater  acoustic 
system  performance,  the  effects  of  internal  waves  on  sound  propagation  in  shallow 
waters  need  to  be  accounted  for. 
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I.  INTRODUCTION 


Internal  waves  have  been  observed  in  many  parts  of  the  world’s  oceans  (Jackson 
2004,  see  Figure  1),  and  sound  waves  propagating  through  the  ocean  are  refracted  by 
these  internal  waves.  In  deep  waters,  such  refraction  generally  only  causes  variability  in 
transmission  time  (Marwoto  2015);  it  does  not  cause  the  sound  waves  to  have  additional 
sea-surface/seabed  boundary  interaction  given  the  depth  of  the  water  column. 


Figure  1.  Locations  of  Observed  Internal  Waves.  Source:  Jackson  (2004). 


The  red  dots  indicate  the  locations  where  internal  waves  are  recorded  in  “An  Atlas  of 
Oceanic  Internal  Solitary  Waves.”  The  lack  of  observed  internal  waves  in  the  open  ocean 
is  due  to  a  lack  of  data  at  those  locations. 


In  shallow  waters  such  as  the  continental  shelves  however,  an  additional 
downward  refracting  of  sound  waves  due  to  internal  waves  could  cause  them  to  interact 
more  often  with  the  seabed,  causing  greater  transmission  loss  than  would  be  expected  in  a 
scenario  without  internal  waves,  as  additional  energy  from  sound  waves  gets  dissipated 
into  or  gets  scattered  by  the  seabed.  Apel  et  al.  (1997,  2007)  have  shown  the  importance 
of  internal  waves  in  shallow  water  acoustics,  while  Zhou  et  al.  (1991)  have  shown  25  dB 
signal  fluctuations  in  the  Yellow  Sea  with  timescales  of  minutes  to  hours. 
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The  objectives  of  this  thesis  are  to  identify  the  types  of  internal  waves  in  the 
shallow  waters  of  a  continental  shelf  and  to  analyze  the  effects  they  have  on  sound 
propagation  ranges.  The  internal  wave  data  were  obtained  from  an  Office  of  Naval 
Research  (ONR)-sponsored  experiment  titled  “2015  Point  Sal  Inner-shelf  Experiment 
(PSIEX),”  which  was  conducted  in  the  coastal  waters  off  Point  Sal,  California.  The 
internal  waves  were  measured  using  vertical  strings  of  temperature  sensors,  while  their 
effects  on  sound  propagation  ranges  were  quantified  via  simulations  using  a  ray  theory 
acoustic  model  called  BELLHOP.  The  end  goal  is  to  do  a  generic  study  to  determine  how 
internal  waves  affect  sound  propagation  ranges  in  the  shallow  waters  of  the  continental 
shelves. 
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II.  PHYSICAL  BACKGROUND 


A.  WHAT  ARE  INTERNAL  WAVES? 

Internal  waves  (IW)  are  gravity  waves  that  oscillate  in  water,  rather  than  on  the 
sea  surface.  An  analogy  would  be  a  layer  of  oil  sitting  on  a  layer  of  water — a  small  shake 
of  the  mixture  will  create  waves  at  the  oil-water  interface.  In  order  for  internal  waves  to 
exist  there  must  be  an  increase  in  sea-water  density  with  increasing  depth,  so  that  any 
perturbation  would  then  cause  an  underwater  wave  to  ripple  outward  as  internal  waves. 

B.  INTERNAL  WAVES  MANIFESTED  AS  WATER  TEMPERATURE 

CHANGES 

In  the  linear  equation  of  state,  density  (p )  is  inversely  proportional  to  temperature 
( T)  and  proportional  to  salinity  (5). 

Linear  Equation  of  State:  — — —  =  -a(T-T0)  +  J3(S  -  S0) , 

Po 

where  p0  =1025  kg/m3,  S0  =  35p  su,  T0=  10  °C,  a  =  2xl0'4  0C4,  /?  =  7.5xl0'4;  a  is  the 
thermal  expansion  coefficient  and  fj  is  the  coefficient  of  saline  contraction. 

For  this  study,  a  constant  salinity  throughout  the  water  column  was  assumed,  as 
there  was  no  introduction  of  fresh  water  into  the  region.  Hence,  temperature  changes 
were  the  only  factor  contributing  to  changes  in  density  in  this  investigation  of  internal 
waves.  As  deeper  waters  are  usually  colder  than  shallower  waters,  this  implies  that  they 
are  also  denser  than  shallower  waters.  Thus,  any  external  forcing  would  set  up  a  wave  at 
the  interface  between  the  waters  of  different  temperature  (and  density).  In  other  words,  an 
internal  wave,  which  propagates  as  perturbations  in  the  density  difference  of  water,  is 
manifested  and  can  be  recorded  as  changes  in  water  temperature. 
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c. 


TYPES  OF  INTERNAL  WAVES  IN  THE  SHALLOW  WATERS  OF  THE 
CONTINENTAL  SHELVES 


1.  Internal  Tides — Internal  Waves  Generated  by  Tidal  Forcing 

Gerkema  and  Zimmerman  (2008)  show  that  internal  waves  on  the  continental 
shelves  are  generated  by  tidal  forcing  over  the  depth  transition  between  the  deep  ocean 
and  the  shallow  continental  shelves.  In  this  scenario,  the  internal  waves  have  a  semi¬ 
diurnal  frequency  of  two  cycles  per  day  (cpd).  This  agrees  well  with  our  experiment 
observations  where  over  a  period  of  one  day,  two  internal  waves  due  to  tidal  forcing  were 
observed  (Figure  2).  Internal  waves  generated  by  tidal  forcing  are  also  called  internal 
tides. 

Very  often,  a  steep  interface  can  occur  between  the  cold  and  warm  front  of  an 
internal  wave.  Such  a  steep  interface  is  called  a  tidal  bore  or,  simply,  bore  (see  Figure  2). 


Figure  2.  Types  of  Internal  Waves  Observed  in  One  Day. 
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2.  Solitons — Solitary  Internal  Waves 

In  addition  to  internal  tides,  shorter  duration  internal  waves  were  observed 
throughout  the  data.  These  short  spikes  of  internal  waves  are  called  internal  solitary 
waves  or  solitons  (Figure  2),  which  are  generated  from  internal  tides  transferring  their 
energies  to  shorter-space  scales  (Gerkema  and  Zimmerman  2008). 
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D.  SEABEDS  AND  THEIR  EFFECTS  ON  SOUND  PROPAGATION 


The  types  of  seabed  in  shallow  waters  play  a  primary  role  in  determining  the 
range  of  sound  propagation,  as  the  waters  in  these  areas  are  usually  downward  refracting 
and  will  cause  sound  to  interact  with  the  seabed  multiple  times  over  short  distances 
(Jensen  et  al.  1994).  Figure  3  illustrates  this  effect. 


Figure  3.  Ray  Trace  of  a  Horizontally  Launched  Ray  in  Mid- water  Column. 


Left:  A  typical  slightly  downward  refracting  sound  speed  profile  recorded  during  the 
experiment  was  used  in  the  ray  trace  simulation  run  as  shown  to  the  right. 


Right:  A  mid-water  column  horizontally  launched  sound  “ray”  gets  refracted  downward 
and  interacts  with  the  seabed  eight  times  in  10  km. 


Three  representative  continental  shelf  seabed  types  were  used  in  the  simulation 
runs:  rocky  seabed  (chalk),  sandy  seabed  (sandy-slit),  and  muddy  seabed  (clay-mud).  The 
detailed  parameters  of  each  bottom  characteristic  used  in  the  BELLHOP  runs  are  listed  in 
Table  1. 


Table  1.  Bottom  Types  Parameters  Used  in  BELLHOP  Runs. 

Adapted  from  Colosi  (2016). 


Bottom  Type 

Compressional 
sound  speed 
(m/s) 

Density 

(kg/m3) 

Attenuation 

(dB/L) 

Rocky  (chalk) 

2400 

2200 

0.2 

Sandy  (sandy-silt) 

1644 

1769 

1.22 

Muddy  (clay-mud) 

1470 

1312 

0.09 
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To  understand  the  effects  of  bottom  interaction  on  acoustics  energy,  it  is 
necessary  to  invoke  the  formula  for  reflection  coefficient  to  determine  how  much  of  the 
original  ray  intensity  was  reflected  from  the  seabed  at  various  grazing  angles. 


m  sin  9  -Jn2-  cos2  9 

Reflection  coefficient:  Rn  = - - — =L  , 

m  sin  0g  +  Jn2  -  cos2  0g 

where  9g  is  the  sound  ray  grazing  angle  as  measured  from  the  seabed  plane, 
P  c 

m  =  — ,  n  =  — ,  p\  and  fn  are  seawater  and  seabed  density,  respectively,  and  c\  and  C2 
P\  C2 

are  sound  speed  in  seawater  and  seabed,  respectively. 

It  should  be  noted  that  while  all  the  seabed  types  are  denser  than  seawater 

'J 

(approximately  1025  kg/m  ),  the  sound  speeds  in  rocky  and  sandy  seabed  are  faster  than 
seawater  (approximately  1500  m/s),  while  the  sound  speed  of  a  muddy  seabed  is  slower 
than  that  of  seawater.  This  has  important  physical  implications — for  the  rocky  and  sandy 
seabed,  total  internal  reflection  (IR12I  =  1)  occurs  when  the  direction  of  sound  propagation 
(or  a  sound  “ray”)  strikes  the  seabed  at  angles,  6g  <  6cr ,  where  6cr  (critical  angle)  is  the 
highest  angle  at  which  total  internal  reflection  can  occur  (Figure  4).  As  for  the  case  of  a 
muddy  seabed,  total  internal  reflection  does  not  occur  as  there  is  no  critical  angle.  Instead 
there  is  what  is  called  a  Brewster  or  intromission  angle  where  |Ri2|  =  0,  in  which  all 
energy  from  the  sound  waves  in  water  gets  dissipated  into  the  seabed  when  0g  is  at  the 
Brewster  angle. 
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Figure  4.  Reflection  Coefficients  of  Rocky,  Sandy,  and  Muddy  Seabeds  for 

the  Various  Grazing  Angles. 


Bottom  Type 

—  Rocky  (chalk) 

— Sandy  (sandy-silt) 
- — Muddy  (clay-mud) 


A  sound  “ray”  is  a  vector  that  is  normal  to  the  sound  wave  front  and  which  traces  the  path 
that  a  sound  wave  travels. 


Referring  to  Figure  4,  the  rocky  seabed  has  the  highest  critical  angle  of 
approximately  52°;  hence,  any  sound  ray  striking  a  rocky  seabed  at  less  than  that  angle 
will  experience  total  internal  reflection.  In  fact,  more  than  55%  of  the  sound  energy  will 
be  reflected  even  when  sound  strikes  the  seabed  at  a  90°  grazing  angle.  This  means  that 
the  rocky  seabed  will  allow  sound  to  travel  the  farthest  as  it  is  the  least  lossy  of  the  three 
seabeds. 

For  a  sandy  seabed,  the  critical  angle  is  much  smaller  compared  to  the  rocky 
seabed  at  approximately  25°.  At  grazing  angles  greater  than  the  critical  angle,  the  portion 
of  reflected  energy  quickly  drops  to  below  40%  at  grazing  angles  of  above  40°.  Hence,  a 
sandy  seabed  is  more  lossy  compared  to  a  rocky  seabed,  and  sound  experiences  more 
attenuation  when  it  interacts  with  a  sandy  seabed. 

For  the  muddy  seabed,  the  percentage  of  energy  reflected  drops  very  rapidly  from 
100%  at  0°  grazing  angle  to  20%  at  6°  grazing  angle.  The  Brewster  angle  is  reached  at 
approximately  12.5°  grazing  angle,  and  all  the  energy  from  the  sound  waves  gets 
absorbed  into  the  seabed.  Beyond  the  Brewster  angle,  up  to  90°  grazing  angle,  only 
approximately  12%  of  the  sound  energy  gets  reflected  back  into  the  water.  Hence,  a 
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muddy  bottom  is  extremely  lossy,  and  sound  will  experience  the  most  attenuation  in  this 
environment. 

The  simulated  transmission  loss  plots  for  the  three  seabed  types  are  as  shown  in 
Figure  5.  The  plots  were  generated  using  the  sound  speed  profile  (SSP)  of  Figure  3,  and 
the  only  difference  between  the  three  plots  is  the  seabed  type.  From  the  transmission  loss 
plots,  it  is  apparent  that  a  rocky  bottom  has  the  least  sound  energy  transmission  loss, 
while  a  sandy  bottom  has  more  transmission  loss  and  a  muddy  bottom  has  the  most 
transmission  loss. 

Figure  5.  Transmission  Loss  Plots  for  a  Rocky  Seabed  (top), 

a  Sandy  Seabed  (center),  and  a  Muddy  Seabed  (bottom). 
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III.  EXPERIMENT  OVERVIEW 


The  data  used  in  this  investigation  was  collected  during  an  ONR-sponsored  effort 
titled  “2015  Point  Sal  Inner-shelf  Experiment  (PSIEX).”  The  objective  of  the  experiment 
was  to  quantify  the  various  oceanographic  phenomena  in  the  coastal  inner-shelf.  The 
experiment  took  place  in  the  waters  off  Point  Sal,  California,  during  the  summer  of  2015. 

The  locations  of  the  various  sensors  are  as  shown  in  Figure  6.  The  experiment 
lasted  52  days  from  15-Jun-2015  07:00  PDT  (yearday  166.29)  to  06-Aug-2015  07:00 
PDT  (yearday  218.29). 


Figure  6.  Experiment  Site  Sensor  Deployment  Locations. 

Source:  Colosi  et  al.  (2016). 
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For  this  study,  the  focus  was  on  temperature  string  050,  as  it  was  the  deepest 
deployed  sensor  of  the  experiment  at  52  m  water  depth.  Temperature  string  050 
consisted  of  five  temperature  sensors  deployed  from  a  depth  of  6.2  m  to  46.1  m  at 
approximately  10  m  spacing,  and  the  water  temperature  was  sampled  once  every  30 
seconds  (or  120  samples  per  hour).  This  allowed  an  internal  wave  frequency  of  up  to  60 
cph  to  be  sampled,  based  on  Nyquist-Shannon  sampling  theorem,  which  states  that  for  a 
waveform  to  be  accurately  sampled,  the  sampling  frequency  must  be  two  times  that  of  the 
sampled  waveform’s  highest  frequency. 

The  maximum  vertical  oscillating  frequency  of  internal  waves  is  limited  by  the 
Brunt-Vaisala  frequency,  or  buoyancy  frequency,  and  any  external  forcing  of  a  water 
column  at  a  frequency  that  is  higher  than  the  Brunt-Vaisala  frequency  will  not  generate 
internal  waves.  Robert  and  William  (1983)  observe  that  the  typical  Brunt-Vaisala 
frequency  in  shallow  water  is  12-15  cph.  Thus,  the  temperature  sensor  sampling  rate  of 
120  samples  per  hour  is  adequate  for  observing  the  internal  waves  expected  in  this 
environment. 
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IV.  CHARACTERISTICS  OF  INTERNAL  WAVES  IN  THE 

EXPERIMENT  AREA 


A.  OVERVIEW  OF  THE  ENTIRE  TEMPERATURE  TIME  SERIES 

The  recordings  of  temperature  string  050  as  a  function  of  depth  and  time  are 
shown  in  Figures  7a-h,  with  Figures  7a-g  each  showing  one  week  of  data  and  Figure  7h 
showing  three  days  of  data  from  week  8. 


Figure  7.  Temperature  Plots  of  Temperature  String  050  for  the 
Entire  Experiment  (Figures  7a-h). 


Figure  7a.  Temperature  Plot — Week  1. 
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Figure  7b.  Temperature  Plot — Week  2. 
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Figure  7c.  Temperature  Plot — Week  3. 
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Figure  7d.  Temperature  Plot — Week  4. 
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Figure  7e.  Temperature  Plot — Week  5. 
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Figure  7f.  Temperature  Plot — Week  6. 
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Figure  7g.  Temperature  Plot — Week  7. 
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Figure  7h.  Temperature  Plot — Week  8,  First  3  Days. 
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B.  INTERNAL  TIDES 

From  these  plots,  it  was  observed  that  there  was  great  variability  in  internal  wave 
activity,  ranging  from  internal  waves  with  small  temperature  difference  (e.g.,  yeardays 
174-177),  to  internal  waves  with  large  temperature  difference  (e.g.,  yeardays  204-208). 
The  internal  waves  were  observed  to  have  two  oscillations  per  day,  matching  the 
semidiurnal  tidal  frequency.  Hence,  it  can  be  concluded  that  internal  waves  in  this  region 
were  largely  due  to  tidal  forcing. 

C.  SOLITONS 

Figure  8  shows  the  number  of  soliton  occurrences  per  day.  Great  variability  in  the 
daily  occurrence  of  solitons  was  also  observed,  ranging  from  no  solitons  during  yeardays 
174  and  178  to  37  solitons  during  yearday  210. 
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Figure  8.  Daily  Occurrence  Rate  of  Solitons.  Source:  Colosi  et  al.  (2016). 


Yearday 


D.  INTERNAL  WAVE  PERIOD,  PROPAGATION  SPEED,  AND  DIRECTION 

Cross-correlation  of  temperature-strings  in  the  experiment  area  (Colosi  et  al. 
2016)  showed  that  the  internal  waves  propagated  west  to  east  at  an  average  speed  of  0.25 
m/s  (Figure  9,  top),  and  at  an  average  bearing  of  80°  with  respect  to  true  north  (Figure  9, 
bottom).  For  ease  of  illustration,  a  purely  eastward  propagation  of  internal  waves  would 
be  assumed  for  the  study,  as  this  is  a  general  analysis  of  how  internal  waves  affect  sound 
propagation  and  is  not  tied  to  a  specific  bathymetry. 
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Figure  9.  Speed  of  Internal  Wave  Propagation  (top)  and  Direction  of  Internal 
Wave  Propagation  (bottom).  Source:  Colosi  et  al.  (2016). 


Internal  tides  in  shallow  waters  are  also  called  tidal  bores,  or  simply  bores.  The  angles 
used  in  the  “Bore  Direction”  plot  are  with  respect  to  east,  and  positive  angles  are  anti¬ 
clockwise  from  east. 


E.  INTERNAL  WAVE  FRONTS  IN  THE  EXPERIMENT  AREA 

The  internal  waves  will  be  assumed  to  have  a  plane  wave  structure  as  they 
propagate  eastward.  This  implies  that  sound  waves  that  travel  eastward  or  westward  will 
cut  across  the  internal  wave  fronts  (red  arrow  of  Figure  10);  hence,  a  range-dependent 
SSP  was  used  to  simulate  the  behaviors  of  these  sound  waves.  As  for  the  sound  waves 
that  travel  northward  or  southward,  they  will  travel  in-parallel  or  along  the  internal  wave 
fronts  (green  arrow  of  Figure  10);  hence,  a  range-independent  sound  speed  profde  was 
used  to  simulate  the  propagation  behavior  of  these  sound  waves. 
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Figure  10.  Sound  Propagation  Along  and  Across  Internal  Wave  Fronts  in  the 

Experiment  Area  (figure  not  drawn  to  scale). 
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Assuming  that  the  internal  waves  propagate  eastward  and  have  a  plane  wave  structure, 
easterly/westerly  sound  propagation  will  cut  across  internal  wave  fronts  (red  arrow)  and 
experience  range-dependent  sound  propagation,  while  northerly/southerly  sound 
propagation  will  run  along  internal  wave  fronts  (green  arrow)  and  experience  range- 
independent  sound  propagation 


F.  TRANSLATING  WATER  TEMPERATURE  TO  SOUND  SPEED 

For  the  translation  of  temperature,  depth,  and  salinity  into  sound  speed,  the  nine- 
term  Mackenzie  equation  was  used  (Mackenzie  1981),  where  sound  speed  (c)  is  a 
function  of  temperature  (T)  in  degrees  Celsius,  salinity  (S)  in  psu,  and  depth  (D)  in 
meters. 
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For  the  nine-term  Mackenzie  equation: 

c(D,S,T)=  1448.96  +  4.591 T  -  5.304  x  lO'2^  +  2.374  x  KT4^  +  1.340  (5-35)  + 
1.630  x  10  2D  +  1.675  x  10  1 D2  -  1.025  x  10  2T(S  -  35)  -  7.139  x  10  nTD3 

the  range  of  the  equation  validity  is  temperature:  2  to  30°C;  salinity:  25  to  40  psu;  depth: 

0  to  8000  m. 

From  the  Mackenzie  equation,  it  can  be  seen  that  in  shallow  waters,  changes  in 
temperature  are  the  primary  cause  of  sound  speed  changes,  while  sound  speed  changes 
due  to  salinity  changes  are  small  in  the  experiment  area  as  there  was  no  introduction  of 
fresh  water  into  the  region.  An  average  salinity  of  34.44  psu  was  used  for  the  simulation 
runs. 

Translation  of  wave  measurements  in  time  into  distance  traveled  could  be 
obtained  by  multiplying  measurement  time  with  internal  wave  propagation  speed  of  0.25 
m/s.  For  semidiurnal  internal  tides  (approximately  12  hours  apart  crest-to-crest),  this 
translates  to  about  10.8  km  between  the  crests  of  the  internal  tides.  Over  two  days,  the 
internal  waves  will  propagate  43.2  km. 

With  the  translation  of  internal  wave  measurements  in  time  to  physical  distance, 
and  the  Mackenzie  equation  to  get  the  SSP,  there  were  sufficient  data  to  generate  a  range- 
dependent  sound  speed  profile  to  analyze  the  effects  of  internal  waves  on  sound 
propagation. 

An  example  of  translating  temperature  readings  into  a  range-dependent  SSP  is 
shown  in  Figure  1 1  for  the  two-day  period  of  yeardays  200-202.  The  temperature  as  a 
function  of  time  and  depth  plot  of  this  particular  period  is  as  shown  in  the  top  plot  of 
Figure  11.  From  this  data,  a  range-dependent  SSP  (with  depth  adjustments  to  facilitate 
BELLHOP  analysis)  was  derived,  as  shown  in  the  bottom  of  Figure  11.  For  the  x-axis, 
the  time  measurements  were  converted  to  range  based  on  the  internal  wave  propagation 
speed,  with  0  km  at  yearday  200,  and  yearday  202  converts  to  43.2  km  west  of 
temperature  string  050  (as  the  internal  wave  fronts  at  that  location  will  reach  temperature 
string  050  two  days  later  with  an  eastward  propagation).  Such  range-dependent  SSPs 
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were  fed  into  BELLHOP  for  simulations  runs,  to  determine  their  effects  on  sound 
propagation. 


Ligure  1 1 .  Temperature  as  a  Function  of  Time  and  Depth  (top) 

Translated  to  a  Range-Dependent  Sound  Speed  Profile  (bottom). 
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It  should  be  noted  that  an  increase  in  temperature  leads  to  an  increase  in  sound 
speed;  hence,  the  structure  of  the  temperature  plot  and  the  range-dependent  SSP  plot  look 
similar.  This  means  that  the  structure  of  the  internal  waves  has  a  direct  influence  on 
sound  speed  and  the  shape  of  the  internal  waves  can  be  seen  from  the  range-dependent 
SSP. 
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V.  BELLHOP— THE  MODELING  AND  SIMULATION  MODEL 

CHOSEN 


A.  RAY  TRACING 

A  ray  theory  acoustic  model  called  BELLHOP  was  used  to  ray  trace  the  sound 
wave  propagation  path  through  the  water.  An  example  of  a  ray  trace  plot  has  been 
previously  illustrated  in  Figure  3.  The  physics  behind  ray  tracing  is  Snell’s  Law,  which 
calculates  the  refraction  of  sound  waves  as  they  move  horizontally  through  waters  of 
different  density  and  hence  sound  speeds.  A  depiction  of  Snell’s  Law  is  shown  in  Figure 
12. 


Figure  12.  Depiction  of  Snell’s  Law.  Adapted  from  Kapolka  (2015). 
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Ci,  C2:  Sound  speed  of  layer  1  and  layer  2,  respectively. 

Xi,X2:  Wavelength  of  sound  wave  in  layer  1  and  layer  2. 

0i,  02 :  Grazing  angle,  angle  of  sound  “ray”  with  respect  to  the  plane  dividing  the  two 
mediums. 


c  c 

Given  that — - —  =  — - — ,  sound  waves  will  refract  downward  as  sound  speed 
cos  6X  cos  02 

decreases  with  depth,  and  will  refract  upward  as  sound  speed  increases  with  depth  (see 
Figure  13). 


19 


Figure  13.  Refraction  of  Sound  Waves.  Adapted  from  Kapolka  (2015). 
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Top:  Sound  wave  refracted  downward  as  sound  speed  decreases  with  depth. 
Bottom:  Sound  wave  refracted  upward  as  sound  speed  increases  with  depth. 
Ci,  C2,  C3,  C4:  Sound  speed  of  the  various  layers. 


To  generate  the  transmission  loss  plots  as  seen  in  Figure  5,  BELLHOP  ray  traces 
a  large  number  of  sound  “rays”  emanating  from  a  sound  source  and  then  sums  the 
intensity  of  the  rays  as  they  move  downrange. 

B.  WHY  BELLHOP  RAY  THEORY  ACOUSTIC  MODEL 

A  ray  theory  acoustic  model  was  chosen  as  this  study  uses  ray  tracing  of  the 
sound  waves  to  demonstrate  how  internal  waves  modulate  their  interactions  with  the 
seabed.  The  other  reason  was  that  a  ray  theory  acoustic  model  could  also  support  range- 
dependent  SSPs  (Etter  2013),  which  are  required  in  this  study  as  internal  waves  will 
cause  the  SSPs  to  vary  with  range. 
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C.  BELLHOP’S  ENVIRONMENTAL  PARAMETERS 

1.  Adjustment  of  Actual  Sea  Depth  to  BELLHOP  Model’s  Depth 

As  there  were  no  temperature  sensors  at  the  sea  surface  and  seabed  for 
temperature  string  050,  a  translational  step  was  made  to  ensure  compatibility  with 
BELLHOP  by  treating  the  shallowest  sensor  as  being  at  the  top  boundary  of  the  water 
column  and  the  deepest  sensor  as  being  at  the  bottom  boundary  of  the  water  column  in 
BELLHOP.  This  eliminates  the  need  to  extrapolate  the  temperature  readings  to  the  sea 
surface  and  the  seabed. 

The  actual  depth  of  the  temperature  sensors  and  their  adjusted  depth  in 
BELLHOP  are  as  listed  in  Table  2.  The  depth  adjustment  was  done  by  treating  the 
shallowest  sensor  depth  (at  6.2  m)  as  the  “surface  sensor”  (of  0  m)  and  adjusting  all  the 
sensor  depths  upward  by  6.2  m,  with  the  depth  of  the  deepest  sensor  rounded  to  40  m. 


Table  2.  Temperature  String  050 — Temperature  Sensor  Depths. 


Actual  sensor  depth 

Adjusted  depth  for 
BELLHOP  run  (by  -6.2 
m) 

6.2  m 

0.0  m 

16.3  m 

10.1  m 

29.7  m 

23.5  m 

36.1  m 

29.9  m 

46.1  m 

Rounded  to  40.0  m 

2.  Flat  Bathymetry  and  Smooth  Sea  Surface 

To  enable  a  general  study  of  the  effect  of  internal  waves  on  sound  propagation,  a 
flat  seabed  and  a  smooth  sea  surface  were  used  in  the  simulation  runs,  so  as  not  to  limit 
the  applicability  of  the  analysis  to  a  specific  location. 
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3. 


Sound  Source  Depth 


The  sound  source  was  placed  at  a  depth  of  38  m  for  the  simulation  runs  so  that 
sound  rays  launched  at  small  angles  with  respect  to  the  horizontal  plane  would  not  be 
refracted  downward  too  much  before  striking  the  seabed.  This  allows  them  to  graze  the 
seabed  at  a  small  angle,  so  that  most  of  their  energy  will  be  reflected,  thus  enabling  them 
to  propagate  farther,  as  is  the  case  for  many  underwater  applications. 

D.  MEASURES  TO  MAXIMIZE  SIMULATION  ACCURACY 

1.  Ray  Trace  Using  a  Frequency  of  2  kHz 

Etter  (2013)  states  that  the  geometrical  acoustics  approximation  made  for  the  ray 
theory  model,  which  assumes  that  sound  speed  does  not  change  much  over  one  wave 
length  horizontally,  limits  the  ray  acoustics  to  higher  frequency  applications  (as  lower 
frequency  sound  waves  have  longer  wave  lengths). 

Etter  (2013)  also  gives  a  guideline  for  defining  the  frequency  above  which  a  ray 
acoustics  model  becomes  suitable  for  use,  in  the  following  equation: 


where  /  is  the  frequency,  H  is  the  water  depth,  and  c  is  the  speed  of  sound.  With  a 
BELLHOP  simulation  depth  of  40  m,  and  using  an  average  sound  speed  of  1500  m/s,  the 
formula  resolves  to  f  >  375  Hz.  A  frequency  of  2  kHz,  which  is  around  the  frequency 
range  of  a  mid-frequency  sonar,  was  used  for  the  study  and  is  well  above  the  375  Hz 
threshold. 


2.  Gaussian  Beam  Option 

Ray  trace  models  are  unable  to  approximate  the  effects  of  caustics,  which  are 
regions  of  intense  beam  focusing.  One  method  to  mitigate  this  is  Gaussian  beam  tracing. 
Porter  and  Bucker  (1987)  state  that  Gaussian  beam  tracing  avoids  the  traditional  ray  trace 
artefact  of  infinite  energy  levels  at  caustics,  and  is  attractive  for  high-frequency,  range- 


22 


dependent  applications.  Hence,  the  Gaussian  beam  option  in  BELLHOP  was  used  to 
increase  the  ray  trace  modeling  accuracy. 

3.  Cubic  Spline  Fitting  of  the  SSP  Plot 

The  depth  interval  between  sensors  on  temperature  string  050  was  10  meters. 
Such  a  coarse  depth  resolution  could  cause  inaccuracies  in  ray  tracing  if  straight  lines 
(linear  interpolation)  were  used  to  join  the  data  points.  To  overcome  this,  the  cubic  spline 
interpolation  of  the  SSP  was  used  in  the  BELLHOP  runs,  which  replaced  the  linear 
approximation  with  a  more  realistic  curved  line  for  the  SSP  (see  Figure  14). 


Figure  14.  Comparison  between  Straight  Line  Linear  Approximation  and 
Cubic  Spline  Interpolation  of  the  Sound  Speed  Profde. 


4.  Incoherent  Beam  Option 

The  BELLHOP  simulation  runs  were  conducted  using  the  “incoherent”  beam 
option,  which  allows  a  smoothed  transmission  loss  plot  of  the  sound  energy  over  range.  A 
“coherent”  beam  option  would  have  accounted  for  the  phase  of  the  sound  waves  and 
shown  the  constructive  and  destructive  interferences  of  the  signal  over  range  as 
distinctive  peaks  and  null  in  the  transmission  loss  plot.  The  difference  in  transmission 
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loss  between  an  incoherent  and  a  coherent  beam  option  is  illustrated  in  Figure  15.  The 
incoherent  beam  option  was  chosen  as  this  study  is  focused  on  the  effects  of  internal 
waves  on  sound  propagation  over  longer  ranges,  and  does  not  require  an  in-depth  analysis 
of  the  effects  of  constructive  and  destructive  interferences  of  the  signal  over  short 
distances  of  a  few  wave  lengths. 


Figure  15.  Incoherent  Transmission  Loss  (top), 
Coherent  Transmission  Loss  (bottom). 


Transmission  loss  plot  with  BELLHOP'S  incoherent  beam  option 


Range  (km) 

Transmission  loss  plot  with  BELLHOP'S  coherent  beam  option 


Transmission  Loss 
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Transmission  Loss 

m 


20  25 
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VI.  TYPES  OF  ANALYSIS  AND  THEIR  CONSIDERATIONS 

A.  TIME  PERIODS  CHOSEN  FOR  ANALYSIS 

Five  time  periods  with  distinct  internal  wave  characteristics  were  chosen  from  the 
experiment  data  for  the  BELLHOP  simulation  runs  (Table  3).  The  first  was  the  period  at 
yeardays  175-177,  which  was  a  period  of  well-defined  internal  tides  with  few  solitons. 
The  difference  between  the  surface  and  bottom  water  temperature  in  this  period  was 
small  at  2.5°C. 

The  second  period  was  at  yeardays  205-207,  which  was  also  a  period  of  well- 
defined  internal  tides,  but  with  many  solitons.  The  surface  and  bottom  water  temperature 
difference  in  this  period  was  large  at  7°C. 

The  third  period  was  at  yeardays  210-212,  which  was  a  period  of  weak  internal 
tides  but  with  a  large  number  of  soliton  occurrences  and  a  surface  and  bottom  water 
temperature  difference  of  6°C. 

The  fourth  period,  at  yeardays  200.6-202.9,  had  well-defined  internal  tides  but 
few  solitons  and  a  surface  and  bottom  water  temperature  difference  of  5°C.  Two 
simulation  runs  were  conducted  for  period  4,  one  with  the  sound  source  in  the  middle  of 
the  cold  front  of  an  internal  tide  (yearday  200.6,  period  4a  of  Table  3),  the  other  with  the 
sound  source  in  the  middle  of  the  warm  front  of  an  internal  tide  (yearday  200.9,  period  4b 
of  Table  3). 

The  fifth  period  was  at  approximately  yearday  174.22,  which  was  at  the  warm  and 
cold  interface  of  a  steep  tidal  bore.  Two  simulation  runs  were  conducted  at  data  points 
that  were  5  minutes  or  75  m  apart,  with  one  run  at  the  warm  front  of  the  tidal  bore  and  the 
other  run  at  the  cold  front  of  the  tidal  bore. 
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Table  3.  Internal  Wave  Characteristics  of  the  Selected  Time  Periods  for  Analysis. 


Period 

No. 


Yearday 


Internal  tide 


Solitons 


Temperature  across  time  and  depth 


1 


175-177 


-Well-defined 

-Small 

temperature 

difference 

(23°C) 


Few 


m  17Bl3  1761  Uii  1TSJ  17ft  176  3  1H4  1TB  6  I7ftft  ttt 


205-201 


-Well-defined 
-Large 
temperature 
difference  (7°C) 


Many 


210-212 


-Weak 
-Large 
temperature 
difference  (6°C) 


Mam 


11(4  ZIDvfl  3101  31 1  3112  3114  2116  211 J  292 


4a 


200.6- 
202.6 
Cold 
front  of 

IW 


-Well-defined 
-Fairly  large 
temperature 
difference  (5°C) 


Few 


l  i 

m 

Ji 

La 

8MB  3»B  JOl  291 J  8011  3SJ1B  Zn* 


mi  to?  i  awfl 


Starts  at  the  cold  front  of  an  internal  tide 


4b 


200.9- 
202.9 
Warm 
front  of 

IW 


Same  as  period 
4b 


Same  as 
period  4b 


3011  2011  301.5  201. T  2019  2021  202  3-  202  5  202  7 


Starts  at  the  warm  front  of  an  internal  tide 


road ay 


Two  runs  at  5  minutes  or  75  m  part  were  conducted; 
one  starts  at  tlie  cold  front  of  the  sleep  tidal  bore,  the 
other  starts  at  the  warm  front  of  the  tidal  bore. 
(Note:  This  plot  has  a  different  color  scale  from  [he 
previous  plots  to  accentuate  the  temperature 
differences) 
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B.  BELLHOP  PLOTS  USED  FOR  ANALYSIS 


Transmission  loss  plots  are  generated  using  range-independent  and  range- 
dependent  SSPs,  so  that  a  quantitative  comparison  could  be  made  on  how  much  internal 
waves  affect  sound  propagation  (range-dependent  SSP)  as  compared  to  a  model  that  does 
not  account  for  internal  waves  (range-independent  SSP).  An  example  of  a  range- 
independent  and  range-dependent  sound  speed  profde  is  shown  in  Figure  16. 


Figure  16.  Range-Independent  SSP  (top),  Range-Dependent  SSP  (bottom). 
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It  has  been  mentioned  that  a  range-independent  SSP  would  be  more  appropriate 
for  simulating  along  internal  wave  front  sound  propagation  (northward/southward  in  the 
experiment  area),  while  the  range-dependent  SSP  would  be  more  appropriate  for 
simulating  across  internal  wave  front  sound  propagation  (eastward/westward  in  the 
experiment  area). 

Two  types  of  transmission  plots  are  generated;  the  first  is  the  “traditional” 
transmission  loss  plot  that  has  been  shown  in  Figure  5,  which  uses  a  gradual  color 
scheme  to  show  the  gradual  loss  of  signal  strength  over  range.  Such  plots  allow  the  sound 
loss  “structure”  to  be  seen. 

The  second  type  of  transmission  loss  plot  is  the  70  dB  transmission  range  plot,  an 
example  of  which  is  shown  in  Figure  17.  Such  plots  allow  a  quantitative  assessment  of 
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the  range  a  sound  signal  can  reach  before  its  transmission  loss  exceeds  70  dB.  Referring 
to  Figure  18,  yellow  shows  the  region  where  transmission  loss  is  less  than  70  dB,  while 
dark  blue  shows  the  region  where  transmission  loss  is  more  than  70  dB.  A  transmission 
loss  of  70  dB  was  chosen  as  the  cut-off  value  in  this  study  as  it  allows  the  display  area  of 
the  plot  to  be  maximized,  which  for  the  case  of  a  rocky  bottom  of  Figure  18,  the  70  dB 
transmission  loss  is  at  40  km,  which  is  almost  to  the  end  of  the  plot.  For  this  study,  the  70 
dB  transmission  ranges  are  measured  from  the  source  range  (0  km)  to  the  maximum 
range  that  the  70  dB  contour  reaches. 


Figure  17.  Example  of  a  70  dB  Transmission  Range  Plot. 

Period  1,  rocky  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 
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Range  (km) 


28 


VII.  RESULTS  OF  MODELING  AND  SIMULATION  RUNS 


A.  PERIOD  1:  INTERNAL  TIDES  WITH  SMALL  TEMPERATURE 

DIFFERENCE,  FEW  SOLITONS 

The  first  period  chosen  for  analysis  (yeardays  175-177)  was  characterized  by 
well-defined  internal  tides,  with  a  small  temperature  difference  of  2.5°C  (from  10  to 
12.5°C)  between  the  bottom  and  the  surface  of  the  water,  and  few  solitons  (fewer  than  5 
solitons  a  day,  see  Figure  8). 

The  sound  speed  profile  ranged  from  approximately  1489  m/s  to  1498  m/s.  The 
range-independent  sound  speed  profile  and  range-dependent  sound  speed  profile  are 
shown  in  Figure  18,  and  the  SSP  at  0  km  is  at  yearday  175.  It  can  be  seen  from  the  range- 
independent  sound  speed  profile  (Figure  19,  top)  that  a  bottom  layer  exists  from  30-40  m 
depth,  and  this  should  serve  to  give  an  extended  transmission  range  in  the  bottom  layer 
due  to  downward  refraction  trapping  sound  energy  in  the  lower  water  column. 

The  transmission  loss  plots  for  the  three  types  of  seabed,  with  a  range- 
independent  SSP  and  range-dependent  SSP,  are  shown  in  Figure  19.  The  70  dB 
transmission  range  plots  are  shown  in  Figure  20.  The  sound  source  depth  for  all  these 
plots  is  38  m. 

Across  all  bottom  types,  the  modeling  runs  that  accounted  for  internal  waves  gave 
a  reduced  transmission  range  compared  to  a  modeling  run  that  does  not  factor  the  effects 
of  internal  waves.  For  the  rocky  seabed,  the  70  dB  transmission  range  dropped  from  40.0 
km  for  along  internal  wave  front  propagation  to  21.5  km  for  across  internal  wave  front 
propagation.  For  the  sandy  seabed,  the  70  dB  transmission  range  dropped  from  27.3  km 
to  6.2  km,  while  for  the  muddy  seabed,  the  70  dB  transmission  range  dropped  from  26.5 
km  to  2.5  km.  Hence,  it  can  be  seen  that  for  a  small  temperature  difference  between  the 
seabed  and  sea  surface,  and  with  a  bottom  layer,  factoring  the  effects  of  internal  waves 
resulted  in  significant  acoustics  range  reduction  compared  to  a  model  that  does  not 
account  for  internal  waves. 
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Figure  18.  Period  1 — Range-Independent  SSP  (top),  Range-Dependent 

SSP  (bottom)  0  km  at  Yearday  175. 
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Figure  19.  Period  1 — 70  dB  Transmission  Range  for  a 

Rocky,  Sandy,  and  Muddy  Seabed.  (Figures  19a-c). 


Figure  19a.  Period  1 — Rocky  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Period  1,  rocky  bottom  transmission  loss  through  a  range-dependent  SSP 
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Figure  19b.  Period  1 — Sandy  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 

Period  1,  sandy  bottom  transmission  loss  through  a  range- in  dependent  SSP  Transmission  loss 
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Figure  19c.  Period  1 — Muddy  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 

Period  1,  muddy  bottom  transmission  loss  through  a  range-independent  SSP  Transmission  Loss 
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Period  1,  muddy  bottom  transmission  loss  through  a  range-dependent  SSP 
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Depth(m)  Depth(m) 


Figure  20.  Period  1 — 70  dB  Transmission  Range  for  a  Rocky,  Sandy,  and 

Muddy  Seabed.  (Figures  20a-c). 


Figure  20a.  Period  1 — Rocky  Bottom  70  dB  Transmission  Range  for  a  Range- 
Independent  SSP  and  a  Range-Dependent  SSP. 

Period  1,  rocky  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 

10 
20 
30 
40 

0  5  10  15  20  25  30  35  40  43.2 

0 
10 
20 
30 
40 

0  5  10  15  20  25  30  35  40  43.2 

Range  (km) 


Range  (km) 

Period  1,  rocky  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front) 
(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 


Figure  20b.  Period  1 — Muddy  Bottom  70  dB  Transmission  Range  for  a  Range- 

Independent  SSP  and  a  Range-Dependent  SSP. 


Period  1,  sandy  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 
(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Period  1.  sandy  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Figure  20c.  Period  1 — Muddy  Bottom  70  dB  Transmission  Range  for  a  Range- 
Independent  SSP  and  a  Range-Dependent  SSP. 

Period  1t  muddy  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 

(Yellow;  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Period  1,  muddy  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front) 
(Yellow;  Tx  loss  <  70  dB;  Blue;  Tx  loss  >  70  dB) 


B.  PERIOD  2:  INTERNAL  TIDES  WITH  LARGE  TEMPERATURE 

DIFFERENCE,  MANY  SOLITONS 

This  period  (yeardays  205-207)  was  characterized  by  well-defined  internal  tides, 
with  a  large  temperature  difference  of  7°C  (from  1 1  to  18°C)  between  the  bottom  and  the 
surface  of  the  water,  and  many  solitons  (more  than  15  solitons  per  day,  see  Figure  8).  The 
sound  speed  profde  ranged  from  approximately  1494  m/s  to  1515  m/s.  The  range- 
independent  sound  speed  profde  and  range-dependent  sound  speed  profde  are  shown  in 
Figure  21,  and  the  SSP  at  0  km  is  centered  at  yearday  205.  The  transmission  loss  plots  for 
the  three  types  of  seabed  are  shown  in  Figure  22,  while  the  70  dB  transmission  range 
plots  are  shown  in  Figure  23. 

Similar  to  period  1,  factoring  for  internal  waves  in  the  model  causes  more 
attenuation  compared  to  a  modeling  run  without  internal  waves.  For  the  rocky  seabed,  the 
70  dB  transmission  range  dropped  from  37.5  km  to  18.3  km  when  the  effects  of  internal 
waves  were  accounted  for.  For  the  sandy  seabed,  the  70  dB  transmission  range  dropped 
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from  12.2  km  to  6.5  km,  while  for  the  muddy  seabed,  the  70  dB  transmission  range 
actually  increased  from  1.8  km  to  2.5  km. 


Figure  21.  Period  2 — Range-Independent  SSP  (top), 

Range-Dependent  SSP  (bottom)  0  km  at  Yearday  205. 
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Figure  22.  Period  2 — Transmission  Loss  for  a  Rocky,  Sandy,  and  Muddy 

Seabed.  (Figures  22a-c). 


Figure  22a.  Period  2 — Rocky  Bottom  Transmission  Loss  through  Range- 

Independent  SSP  and  Range-Dependent  SSP. 
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Figure  22b.  Period  2 — Sandy  Bottom  Transmission  Loss  through 

Range-Independent  SSP  and  Range-Dependent  SSP. 
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Period  2,  sandy  bottom  transmission  loss  through  a  range-dependent  SSP 
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Figure  22c.  Period  2 — Muddy  Bottom  Transmission  Loss  through 

Range-Independent  SSP  and  Range-Dependent  SSP. 
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Figure  23.  Period  2 — 70  dB  Transmission  Range  for  a  Rocky,  Sandy,  and 

Muddy  Seabed.  (Figures  23a-c). 

Figure  23  a.  Period  2 — Rocky  Bottom  70  dB  Transmission  Range  for  a 

Range-Independent  SSP  and  a  Range-Dependent  SSP. 

Period  2,  rocky  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Period  2,  rocky  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front) 
(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 


Figure  23b.  Period  2 — Sandy  Bottom  70  dB  Transmission  Range  for  a 

Range-Independent  SSP  and  a  Range-Dependent  SSP. 


Period  2,  sandy  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front) 
(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Figure  23c.  Period  2 — Muddy  Bottom  70  dB  Transmission  Range  for  a 

Range-Independent  SSP  and  a  Range-Dependent  SSP. 

Period  2,  muddy  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 
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Period  2,  muddy  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front} 


Range  (km) 


C.  PERIOD  3:  WEAK  INTERNAL  TIDES,  MANY  SOLITONS 

This  period  (yeardays  210-212)  was  characterized  by  weak  internal  tides  and 
many  solitons  (more  than  15  solitons  per  day,  see  Figure  8).  The  water  temperature 
difference  was  strong  at  6°C  (from  1 1.5  to  17.5°C)  between  the  bottom  and  the  surface  of 
the  water.  The  sound  speed  profile  ranged  from  approximately  1495  m/s  to  1514  m/s.  The 
range-independent  sound  speed  profile  and  range-dependent  sound  speed  profile  are 
shown  in  Figure  24,  and  the  SSP  at  0  km  is  at  yearday  210.  The  transmission  loss  plots 
for  the  three  types  of  seabed  are  shown  in  Figure  25,  while  the  70  dB  transmission  range 
plots  are  shown  in  Figure  26. 

The  simulation  showed  that  even  with  internal  waves  only  in  the  form  of  solitons, 
they  still  caused  a  range  reduction  as  compared  to  a  simulation  that  does  not  account  for 
internal  waves.  For  a  rocky  bottom,  the  70  dB  transmission  range  was  reduced  from  36.8 
km  for  along  internal  wave  front  propagation  to  18.5  km  for  across  internal  wave  front 
propagation.  For  a  sandy  bottom,  the  range  was  reduced  from  17.2  km  to  9.5  km.  For  a 
muddy  bottom,  the  range  was  reduced  from  7.5  km  to  3.7  km. 
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Figure  24.  Period  3 — Range-Independent  SSP  (top), 

Range-Dependent  SSP  (bottom)  0  km  at  Yearday  210. 
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Figure  25.  Period  3 — Transmission  Loss  for  a  Rocky,  Sandy,  and  Muddy 

Seabed.  (Figures  25a-c). 


Figure  25a.  Period  3 — Rocky  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Figure  25b.  Period  3 — Sandy  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Figure  25c.  Period  3 — Muddy  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Period  3,  muddy  bottom  transmission  loss  through  a  range-dependent  SSP 


0  5  10  15  20  25  30  35  40  43,2 

Range  (km) 


39 


Figure  26.  Period  3 — 70  dB  Transmission  Range  for  a  Rocky,  Sandy,  and 

Muddy  Seabed.  (Figures  26a-c). 


Figure  26a.  Period  3 — Rocky  Bottom  70  dB  Transmission  Range  for  a  Range- 
Independent  SSP  and  a  Range-Dependent  SSP. 

Period  3,  rooky  bottom  TOdB  transmission  range,  range -in dependent  SSP  [along  IW  front} 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Figure  26b.  Period  3 — Sandy  Bottom  70  dB  Transmission  Range  for  a  Range- 
Independent  SSP  and  a  Range-Dependent  SSP. 

Period  3,  sandy  bottom  TOdB  transmission  range,  range -independent  SSP  (along  I W  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Figure  26c.  Period  3 — Muddy  Bottom  70  dB  Transmission  Range  for  a 

Range-Independent  SSP  and  a  Range-Dependent  SSP. 

Period  3,  muddy  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 
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Period  3,  muddy  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front} 
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D.  PERIOD  4:  INTERNAL  TIDES  WITH  LARGE  TEMPERATURE 

DIFFERENCE,  FEW  SOLITONS 

Period  4  (yeardays  200.6-202.9)  was  characterized  by  well-defined  internal  tides, 
with  a  relatively  large  temperature  difference  of  5.5°C  (from  10.5  to  16°C)  between  the 
bottom  and  the  surface  of  the  water,  and  few  solitons  (fewer  than  5  solitons  a  day,  refer  to 
Figure  8).  The  sound  speed  profile  ranged  from  approximately  1492  m/s  to  1508  m/s.  For 
period  4a  (yeardays  200.6-202.6),  the  sound  source  was  placed  at  the  center  of  a  cold 
front  of  an  internal  tide  (Figure  27),  while  period  4b  (yeardays  200.9-202.9)  had  the 
sound  source  placed  in  the  center  of  the  warm  front  of  the  same  internal  tide  (Figure  30). 
The  objective  was  to  determine  if  there  are  significant  differences  in  sound  attenuation 
when  a  sound  source  is  in  the  different  temperature  fronts  of  an  internal  tide. 

For  period  4a,  the  transmission  loss  plots  for  the  three  types  of  seabed  are  shown 
in  Figure  28.  The  70  dB  transmission  range  plots  are  shown  in  Figure  29.  For  a  rocky 
bottom,  the  70  dB  transmission  range  was  reduced  from  39.2  km  for  along  internal  wave 
front  propagation  to  18.3  km  for  across  internal  wave  front  propagation.  For  a  sandy 
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bottom,  the  range  was  reduced  from  19.5  km  to  6.7  km.  For  a  muddy  bottom,  the  range 
was  reduced  from  4.3  km  to  2.0  km. 

For  period  4b,  the  transmission  loss  plots  for  the  three  types  of  seabed  are  shown 
in  Figure  31.  The  70  dB  transmission  range  plots  are  shown  in  Figure  32.  For  a  rocky 
bottom,  the  70  dB  transmission  range  was  reduced  from  37.8  km  for  along  internal  wave 
front  propagation  to  19.2  km  for  across  internal  wave  front  propagation.  For  a  sandy 
bottom,  the  range  was  reduced  from  10.0  km  to  5.8  km.  For  a  muddy  bottom,  the  range 
remains  the  same  at  2.0  km. 

It  can  be  observed  that  for  a  sandy  and  a  muddy  seabed,  having  a  sound  source  in 
the  cold  or  warm  front  of  an  internal  tide  results  in  a  large  variability  in  along  internal 
wave  front  sound  propagation  ranges.  For  a  sandy  bottom,  the  along  cold  front  70  dB 
transmission  range  was  19.5  km,  while  the  along  warm  front  70  dB  transmission  range 
was  10.0  km.  For  a  muddy  bottom,  the  along  cold  front  70  dB  transmission  range  was  4.3 
km,  while  the  along  warm  front  70  dB  transmission  range  was  2.0  km. 


Figure  27.  Period  4a — Range-Independent  SSP  (top),  Range-Dependent  SSP 
(bottom)  0  km  at  Yearday  200.6  (Cold  Front  of  IW). 
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Figure  28.  Period  4a — Transmission  Loss  for  a  Rocky,  Sandy,  and  Muddy 

Seabed.  (Figures  28a-c). 


Figure  28a.  Period  4a — Rocky  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 

Period  4a,  rocky  bottom  transmission  loss  through  a  range-independent  SSP  Transmission  Loss 


0  5  10  15  20  25  30  35  40  43.2 

Range  (km) 

Period  4a,  rocky  bottom  transmission  loss  through  a  range-dependent  SSP 


0  5  10  15  20  25  30  35  40  43.2 

Range  (km) 


Figure  28b.  Period  4a — Sandy  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Figure  28c.  Period  4a — Muddy  Bottom  Transmission  Loss  through  Range- 

Independent  SSP  and  Range-Dependent  SSP. 


Period  4a,  muddy  bottom  transmission  loss  through  a  range-independent  SSP  Transmission  Loss 


Range  (km) 


Period  4a,  muddy  bottom  transmission  loss  through  a  range-dependent  SSP 
Across  internal  wave  front  sound  propagation 


0  5  10  15  20  25  30  35  40  43.2 


I 


20 

40 


60 


I 


80 


Range  (km) 


Figure  29.  Period  4a — 70  dB  Transmission  Loss  Range  for  a  Rocky,  Sandy,  and 

Muddy  Seabed.  (Figures  29a-c). 


Figure  29a.  Period  4a — Rocky  Bottom  70  dB  Transmission  Range  for  a  Range - 

Independent  SSP  and  a  Range-Dependent  SSP. 

Period  4af  rocky  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 
_  (Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Range  (km) 


432 


43.2 


44 


Depth(m)  Depth(m)  Depth(m)  Depth(m) 


Figure  29b.  Period  4a — Sandy  Bottom  70  dB  Transmission  Range  for  a  Range- 

Independent  SSP  and  a  Range-Dependent  SSP. 


Period  4a,  sandy  bottom  70dB  transmission  range,  range-independent  SSP  (along  1W front) 
(Yellow;  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Period  4a,  sandy  bottom  70dB  transmission  range,  range-dependent  SSP  (across  IW  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue;  Tx  loss  >  70  dB) 
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Figure  29c.  Period  4a — Muddy  Bottom  70  dB  Transmission  Range  for  a  Range- 

Independent  SSP  and  a  Range-Dependent  SSP. 

Period  4a,  muddy  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Figure  30.  Period  4b — Range-Independent  SSP  (top),  Range-Dependent  SSP 
(bottom)  0  km  at  Yearday  200.9  (Warm  Front  of  IW). 
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Figure  3 1 .  Period  4b — Transmission  Loss  for  a  Rocky,  Sandy,  and  Muddy 

Seabed.  (Figures  31a-c). 


Figure  31a.  Period  4b — Rocky  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Figure  3  lb.  Period  4b — Sandy  Bottom  Transmission  Loss  through  Range-Independent 

SSP  and  Range-Dependent  SSP. 
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Figure  31c.  Period  4b — Muddy  Bottom  Transmission  Loss  through  Range- 

Independent  SSP  and  Range-Dependent  SSP. 
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Figure  32.  Period  4b — Transmission  Loss  Range  for  a  Rocky,  Sandy,  and 

Muddy  Seabed.  (Figures  32a-c). 


Figure  32a.  Period  4b — Rocky  Bottom  70  dB  Transmission  Range  for  a  Range- 

Independent  SSP  and  a  Range-Dependent  SSP. 

Period  4b,  rocky  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 
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Period  4b — Sandy  Bottom  70  dB  Transmission  Range  for  a  Range- 
Independent  SSP  and  a  Range-Dependent  SSP. 


Period  4b,  sandy  bottom  70dB  transmission  range,  range-independent  SSP  (along  IW  front) 

(Yellow:  Tx  loss  <  70  dB;  Blue:  Tx  loss  >  70  dB) 
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Figure  32c.  Period  4b — Muddy  Bottom  70  dB  Transmission  Range  for  a  Range- 

Independent  SSP  and  a  Range-Dependent  SSP. 
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E.  PERIOD  5:  IN  THE  WARM  AND  COLD  FRONT  OF  A  STEEP  TIDAL 

BORE 

Period  5  was  at  approximately  yearday  174.22,  at  the  warm  and  cold  interface  of  a 
steep-face  tidal  bore.  Two  simulation  runs  were  conducted  at  data  points  that  were  5 
minutes  or  75  m  apart,  with  one  data  point  at  the  warm  front  of  the  tidal  bore  and  the 
other  at  the  cold  front  of  the  tidal  bore.  The  exact  time  of  analysis  and  the  70  dB 
transmission  range  for  a  range-independent  SSP  and  a  range-dependent  SSP  are  shown  in 
Figure  33  and  Figure  34,  with  Figure  31  starting  in  the  warm  front  and  Figure  32  starting 
in  the  cold  front.  A  sandy  bottom  was  used  in  the  analysis. 

While  the  range-dependent  SSP  (across  internal  wave  fronts)  did  not  show  much 
change  in  the  70  dB  transmission  range  between  the  two  simulation  time  periods  (8.3  km 
and  8.9  km),  it  was  a  different  story  for  the  range-independent  SSP  (along  internal  wave 
fronts).  For  the  range-independent  SSP  the  coming  of  the  cold  front  of  the  internal  tidal 
bore  to  displace  the  warmer  waters  resulted  in  the  70  dB  transmission  range  increasing 

from  12.1  km  to  32.1  km  in  just  5  minutes  or  just  75  m  apart. 
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Figure  33.  Period  5 — Sandy  Bottom  70  dB  Transmission  Range  at  the  Warm 

Front  of  a  Tidal  Bore. 
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Figure  34.  Period  5 — Sandy  Bottom  70  dB  Transmission  Range  at  the  Cold 

Front  of  a  Tidal  Bore. 


Internal  wave  profile,  0  km  range  at  yearday  174.2292  (23-Jur-201S,  Q5:30h) 
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F.  SUMMARY  AND  INTERPRETATION  OF  SIMULATION  RESULTS 

The  results  of  the  transmission  loss  runs  are  summarized  in  Table  4. 


Table  4.  Transmission  Ranges  for  the  Various  Time  Periods. 


Period 

Centered 

yearday 

Internal  tide 

Solitons 

70  dB  transmission  loss  range 

1 

175-177 

-Well-defined 

Few 

Seabed 

Along  IW 

Across  IW 

%  change 

-Small 

Rocky 

40.0  km 

21.5  km 

-46% 

temperature 

Sandy 

27.3  km 

6.2  km 

-77% 

difference 

Muddy 

26.5  km 

2.5  km 

-91% 

(2.5°C) 

2 

205-207 

-Well-defined 

Many 

Seabed 

Along  IW 

Across  IW 

%  change 

-Large 

Rocky 

37.5  km 

18.3  km 

-51% 

temperature 

Sandy 

12.2  km 

6.5  km 

-47% 

difference 

Muddy 

1.8  km 

2.5  km 

39% 

(7°C) 

3 

210-212 

-Weak 

Many 

Seabed 

Along  IW 

Across  IW 

%  change 

-Large 

Rocky 

36.8  km 

18.5  km 

-50% 

temperature 

Sandy 

17.2  km 

9.5  km 

-45% 

difference 

Muddy 

7.5  km 

3.7  km 

-51% 

(6°C) 

4a 

200.6- 

-Well-defined 

Few 

Seabed 

Along  IW 

Across  IW 

%  change 

202.6 

-Fairly  large 

Rocky 

39.2  km 

18.3  km 

-53% 

Cold  front 

temperature 

Sandy 

19.5  km 

6.7  km 

-66% 

of  IW 

difference 

Muddy 

4.3  km 

2.0  km 

-53% 

(5°C) 

4b 

200.9- 

Same  as  period 

Same  as 

Seabed 

Along  IW 

Across  IW 

%  change 

202.9 

4b 

period  4b 

Rocky 

37.8  km 

19.2  km 

-49% 

Warm 

Sandy 

10.0  km 

5.8  km 

-42% 

front  of 

Muddy 

2.0  km 

2.0  km 

0% 

IW 

5 

174.2- 

Same  as  period 

Same  as 

Seabed 

Along  IW 

Across  IW 

%  change 

176.2 
Edge  of 

1 

period  1 

Sandy 

(warm) 

12.0  km 

8.3  km 

-31% 

tidal  bore 

Sandy 

(cold) 

32.0  km 

8.5  km 

-73% 

Referring  to  Table  5,  for  a  rocky  bottom,  sound  traveling  across  internal  wave 
fronts  (range-dependent  SSP)  had  its  70  dB  transmission  loss  range  reduced  between 
46%  to  53%,  when  compared  to  sound  traveling  along  the  internal  wave  front  (range- 
independent  SSP). 
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The  along  internal  wave  70  dB  transmission  ranges  were  from  36.8  km  to  40.0  km 
across  the  various  time  periods,  giving  a  variability  of  3.2  km.  The  across  internal  wave 
70  dB  transmission  ranges  were  from  18.3  km  to  21.5  km  across  the  various  time  periods, 
giving  a  variability  of  3.2  km. 


Table  5.  70  dB  Transmission  Ranges  for  a  Rocky  Bottom. 


Period 

Unique  characteristics 

Along  IW 

Across  IW 

%  change 

1 

-Sound  source  in  a  well-defined  bottom  cold  layer 
-Internal  tides  only 

40.0  km 

21.5  km 

-46% 

2 

-Internal  tides  and  solitons 

37.5  km 

18.3  km 

-51% 

3 

-Solitons  only 

36.8  km 

18.5  km 

-50% 

4a 

-Sound  source  in  the  cold  front  of  an  internal  tide 

39.2  km 

18.3  km 

-53% 

4b 

-Sound  source  in  the  warm  front  of  an  internal  tide 

37.8  km 

19.2  km 

-49% 

Variability: 

3.2  km 

3.2  km 

Referring  to  Table  6,  for  a  sandy  bottom,  sound  traveling  across  internal  wave 
fronts  (range-dependent  SSP)  had  its  70  dB  transmission  loss  range  reduced  between 
31%  to  77%,  when  compared  to  sound  traveling  along  the  internal  wave  front  (range- 
independent  SSP). 

The  along  wave  front  70  dB  transmission  ranges  were  from  10.0  km  to  32.0  km, 
giving  a  variability  of  22  km,  while  the  across  wave  front  transmission  ranges  were  from 
5.8  km  to  9.5  km,  giving  a  variability  of  3.7  km.  It  can  be  observed  that  a  sound  source 
placed  in  the  cold  front  of  an  internal  wave  supported  an  extended  transmission  range  as 
compared  to  a  sound  source  placed  in  the  warm  front  of  an  internal  wave,  and  this 
resulted  in  a  large  variability  in  the  along  internal  wave  front  70  dB  transmission  range 


Table  6.  70  dB  Transmission  Ranges  for  a  Sandy  Bottom. 


Period 

Unique  characteristics 

Along  IW 

Across  IW 

%  change 

1 

-Sound  source  in  a  well-defined  bottom  cold  layer 
-Internal  tides  only 

27.3  km 

6.2  km 

-77% 

2 

-Internal  tides  and  solitons 

12.2  km 

6.5  km 

-47% 

3 

-Solitons  only 

17.2  km 

9.5  km 

-45% 

4a 

-Sound  source  in  the  cold  front  of  an  internal  tide 

19.5  km 

6.7  km 

-66% 

4b 

-Sound  source  in  the  warm  front  of  an  internal  tide 

10.0  km 

5.8  km 

-42% 

5a 

-Sound  source  in  the  warm  front  of  a  tidal  bore 

12.1  km 

8.3  km 

-31% 

5b 

-Sound  source  in  the  cold  front  of  a  tidal  bore 

32.0  km 

8.5  km 

-73% 

Variability: 

22.0  km 

3.7  km 
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Referring  to  Table  7,  for  a  muddy  bottom,  the  along  internal  wave  front  70  dB 
transmission  ranges  were  from  2.0  km  to  26.5  km  (24.5  km  variability),  while  the  across 
wave  front  transmission  ranges  were  from  2.0  km  to  3.7  km  (1.7  km  variability).  The 
only  instance  where  the  muddy  seabed  could  support  a  longer  range  70  dB  transmission 
range  was  in  period  1 ,  which  had  a  thin  well-defined  cold  bottom  layer.  If  this  period  was 
removed  from  the  analysis,  the  muddy  bottom  along  wave  front  70  dB  transmission 
ranges  would  be  from  2.0  km  to  7.5  km  (5.5  km  variability). 

As  for  percentage  change  in  the  70  dB  transmission  range  for  the  across  internal 
wave  fronts  (range-dependent  SSP)  sound  propagation  as  compared  to  the  along  internal 
wave  fronts  (range-independent  SSP),  a  range  gain  of  39%  (from  1.8  km  to  2.5  km)  was 
observed  in  period  2.  This  range  gain  only  occurred  when  the  along  wave  front  sound 
propagation  ranges  were  short  (at  approximately  2  km).  The  range  gain  could  occur  as 
when  the  grazing  sound  wave  angle  increases  from  beyond  the  Brewster  angle,  where  all 
sound  energy  is  absorbed  by  the  muddy  seabed,  there  will  be  some  energy  that  will  be 
reflected  back  into  the  water.  All  the  other  cases  had  range  reduction  from  0-91%. 


Table  7.  70  dB  Transmission  Ranges  for  a  Muddy  Bottom. 


Period 

Unique  characteristics 

Along  IW 

Across  IW 

%  change 

1 

-Sound  source  in  a  well-defined  thin  bottom  cold  layer 
-Internal  tides  only 

26.5  km 

2.5  km 

-91% 

2 

-Internal  tides  and  solitons 

1.8  km 

2.5  km 

39% 

3 

-Solitons  only 

7.5  km 

3.7  km 

-51% 

4a 

-Sound  source  in  the  cold  front  of  an  internal  tide 

4.3  km 

2.0  km 

-53% 

4b 

-Sound  source  in  the  warm  front  of  an  internal  tide 

2.0  km 

2.0  km 

0% 

Variability: 

24.5  km 
(5.5  km  if 
period  1  is 
excluded) 

1.7  km 
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VIII.  HOW  INTERNAL  WAVES  AFFECT  SOUND  TRAVEL  PATHS 


Observing  that  simulation  runs  that  factored  in  the  effects  of  internal  waves 
generally  resulted  in  greater  transmission  loss  compared  to  simulation  runs  that  do  not 
account  for  the  effects  of  internal  waves,  it  would  be  logical  to  conclude  that  internal 
waves  cause  the  sound  waves  to  refract  downward  more.  This  leads  them  to  strike  the 
seabed  at  a  larger  grazing  angle,  resulting  in  fewer  “rays”  experiencing  total  internal 
reflection;  for  “rays”  beyond  the  critical  angle,  it  pushes  them  into  a  steeper  grazing  angle 
with  a  resulting  smaller  reflection  coefficient. 

To  verify  this  hypothesis,  a  ray  trace  was  done  for  period  4a,  at  a  launch  angle  of 
1°  below  the  horizontal  plane,  as  shown  in  Figure  35.  The  top  ray  trace  was  without  the 
effects  of  internal  waves,  and  the  ray  in  the  positive  ranges  (0  to  20  km)  had  15  bottom 
bounces  in  the  positive  ranges  before  it  exited  the  plot.  For  the  lower  ray  trace  plot,  which 
accounted  for  internal  waves,  the  ray  in  the  positive  ranges  had  17  bottom  bounces  before 
it  exited  the  plot,  which  is  two  more  bottom  interactions  than  the  ray  trace  without 
internal  waves.  Comparing  the  internal  wave  ray  trace  plot  with  the  SSP  plot  (Figure  35, 
bottom),  it  can  be  seen  that  when  the  sound  wave  passes  through  the  warm  front  of  an 
internal  wave,  at  the  5-10  km  region,  with  an  associated  larger  temperature  difference, 
sound  waves  get  “depressed”  more  and  strike  the  seabed  at  a  faster  rate  than  when 
passing  through  a  colder  mass  of  water.  This  would  cause  more  transmission  loss  if  the 
grazing  angle  is  greater  than  the  critical  angle  for  which  total  internal  reflection  can  occur 
for  a  rocky  and  sandy  bottom. 

As  for  the  ray  in  the  negative  ranges  of  the  plot  (0  to  -20  km),  the  top  ray  trace 
had  15  bounces,  while  the  bottom  ray  trace  had  14  bottom  bounces,  which  was  one  fewer 
bottom  bounces.  Internal  waves  still  caused  more  losses  in  this  scenario  as  the  ray  trace 
without  internal  waves  had  the  ray  refracted  downward  at  a  constant  depth  of  20  m,  while 
for  the  ray  trace  with  internal  waves,  the  ray  got  refracted  downward  mostly  at  shallower 
depths  of  less  than  20  m.  This  means  that  the  ray  that  experiences  internal  waves  would 
have  had  a  greater  grazing  angle  when  striking  the  bottom,  which  translates  to  more 
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losses  overall,  even  when  it  had  one  fewer  bottom  interactions  compared  to  the  ray  trace 
without  the  effects  of  internal  waves. 


Figure  35.  Ray  Trace  without  Internal  Waves  (top),  Ray  Trace  with  Internal 
Waves  (middle),  and  Temperature  Profile  (bottom). 
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In  conclusion,  ray  tracing  of  sound  waves  confirmed  the  hypothesis  that  internal 
waves  cause  sound  waves  to  refract  downward  more,  causing  more  bottom  interactions 
over  the  same  distance  or  increasing  the  angle  of  grazing  when  striking  the  seabed.  The 
net  effect  was  more  losses  over  the  same  distance  as  compared  to  a  model  without 
internal  waves.  This  would  explain  why  simulations  that  accounted  for  the  effects  of 
internal  waves  via  range-dependent  SSPs  resulted  in  higher  transmission  loss  as 
compared  to  simulations  that  did  not  account  for  the  effects  of  internal  waves. 
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IX.  CONCLUSION 


The  simulation  of  the  effects  of  internal  waves  on  sound  propagation,  based  on 
real  world  internal  wave  data,  showed  that  internal  waves  introduced  great  variability  in 
sound  propagation  ranges  in  the  shallow  waters  of  the  continental  shelves,  as  the  internal 
waves  modulate  the  SSP  in  time  scales  as  short  as  a  few  minutes. 

Assuming  internal  waves  in  the  experiment  area,  which  are  propagating  in  the 
west-east  direction,  have  plane  wave  structures,  sound  waves  propagating  in  a 
northerly/southerly  direction  will  then  run  along  or  in  parallel  to  the  wave  fronts  of  these 
internal  waves  (see  Figure  36).  In  this  scenario,  the  location  of  the  sound  source  with 
respect  to  the  warm  or  cold  front  of  an  internal  wave  front  has  a  great  effect  on  the  sound 
transmission  range,  with  a  cold  front  allowing  a  farther  transmission  range  compared  to  a 
warm  front.  In  one  of  the  simulation  runs  where  the  warm  and  cold  front  interface  of  an 
internal  wave  was  in  the  form  of  a  steep  tidal  bore,  a  transmission  range  of  32.0  km  was 
reached  in  the  cold  front  of  the  bore,  but  dropped  to  12.0  km  in  the  warm  front  of  the 
bore,  which  was  in  a  data  point  just  5  minutes  or  75  m  apart. 
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Figure  36.  Along  Wave  Front  Variability  in  Sound  Transmission 

Range  Due  to  the  Warm  and  Cold  Fronts  of  Internal  Waves. 
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Sound  waves  propagating  along  the  warm  front  of  an  internal  wave  have  greater 
attenuation  compared  to  sound  waves  propagating  along  the  cold  front  of  an  internal 
wave. 


For  sound  waves  propagating  in  an  easterly/westerly  direction,  they  will  cut 
across  the  wave  fronts  of  internal  waves,  resulting  in  them  having  shorter  transmission 
ranges  as  compared  to  sound  waves  of  the  same  power  level  propagating  in  the 
northerly/southerly  direction  (i.e.,  in  parallel  to  the  wave  fronts  of  the  internal  waves). 
This  scenario  is  illustrated  in  Figure  37.  In  simulation  runs  conducted  using  a  sandy 
seabed,  this  resulted  in  a  31-77%  reduction  in  transmission  ranges  (from  10.0-32.0  km  to 
5. 8-9. 5  km)  when  compared  to  sound  waves  that  travel  along  the  internal  wave  fronts. 
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Figure  37.  Across  Wave  Front  Sound  Transmission 
Range  Reduction  Effects  of  Internal  Wave. 
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Sound  waves  have  greater  attenuation  when  propagating  across  internal  wave  fronts 
when  compared  to  propagating  along  internal  wave  fronts. 


In  conclusion,  for  sound  propagation  along  internal  wave  fronts,  internal  waves 
cause  great  variability  in  propagation  ranges,  with  cold  fronts  allowing  farther 
propagation  range  than  warm  fronts.  As  for  sound  propagation  across  internal  wave 
fronts,  internal  waves  cause  them  to  have  a  shorter  propagation  range  compared  to  the 
sound  waves  of  the  same  power  level  that  are  traveling  along  an  internal  wave  front. 

Hence,  the  current  practice  of  dropping  an  XBT  or  CDT  into  one  spot  to  get  the 
SSP  for  that  location,  and  using  that  one  SSP  for  “range  of  the  day”  estimation  of 
acoustical  systems  performance,  with  the  underlying  assumption  that  the  SSP  does  not 
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change  with  location,  range,  or  time  may  not  be  the  best  approach  in  the  shallow  waters 
of  the  continental  shelves. 

Going  forward,  in  scenarios  where  an  accurate  estimation  of  acoustical  system 
range  performance  is  important  for  “mission  success,”  such  as  in  the  search  for  and 
location  of  an  aircraft  “black  box”  that  is  suspected  to  be  on  the  seabed,  the  effects  of 
internal  waves  on  the  “black  box”  underwater  sound  beacon  transmission  range  would 
have  to  be  accounted  for  in  the  modeling  and  simulation  runs  during  “mission  planning,” 
to  ensure  there  are  no  inadvertent  gaps  in  the  search  operation. 

This  study  conducted  analyses  of  the  effects  of  internal  waves  on  sound 
propagation  using  the  assumption  that  internal  wave  fronts  take  on  the  structure  of  plane 
waves.  While  this  allowed  a  baseline  generic  assessment  of  the  effects  of  internal  waves 
on  sound  propagation,  internal  waves  in  the  shallow  waters  of  the  continental  shelves 
have  been  observed  to  have  more  complex  structures  that  are  closely  tied  to  an  area’s 
bathymetry  (Jackson  2004).  Hence,  to  accurately  account  for  the  effects  of  internal  waves 
on  sound  propagation  ranges,  further  data  measurements  are  required  at  the  locations 
where  underwater  acoustical  operations  have  to  be  carried  out. 
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